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Exact characterization of 0(n) tricriticality in two dimensions 
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We propose exact expressions for the conformal anomaly and for three critical exponents of the 
tricritical 0(n) loop model as a function of n in the range —2 < n < 3/2. These findings are based 
on an analogy with known relations between Potts and 0(n) models, and on an exact solution of 
a 'tri-tricritical' Potts model described in the literature. We verify the exact expressions for the 
tricritical 0(n) model by means of a finite-size scaling analysis based on numerical transfer-matrix 
calculations. 
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While exact results exist for a rich collection of uni- 
versality classes of two-dimensional phase transitions, in- 
cluding g-state Potts criticality and tricriticality, and 
0{n) criticality, such results are still absent for the tricrit- 
ical 0(n) model, except for isolated points at n = and 
1. The purpose of the present work is to fill in this gap 
and to provide exact formulas for the conformal anomaly 
and the main scaling dimensions of the tricritical 0{n) 
model as a function of n. These results are not rigorous 
in the mathematical sense but they may still be assumed 
to be exactly true, as we shall argue below. 

The 0(n) model is defined in terms of n-component 
spins on a lattice, with an isotropic pair coupling of the 
form Eij = e[Si ■ Sj) where i and j denote two neighbor- 
ing lattice sites and e is a function. This model can be 
represented by a graph expansion Q , in which n assumes 
the role of a continuous parameter. For this purpose, it is 
especially useful to choose the model on the honeycomb 
lattice, and the function e as e{p) = — log(l -I- xp), where 
a; is a measure of the inverse temperature. Then the 
graph expansion reduces to a gas of non-intersecting and 
non-overlapping loops on the honeycomb lattice . This 
loop gas representation enables further mappings on the 
Kagome 6- vertex model and the Coulomb gas, and has 
therefore played a crucial role as a step in the derivation 
of exact results for the honeycomb 0{n) model 

Just as for the Potts model, tricriticality can be in- 
duced by introducing a sufficient number of vacancies. 
This was already known Q for the case n = which 
describes the collapse of a polymer at the so-called theta 
point, induced by attractive interactions between the 
polymer segments, and for the Ising case ?^ = 1 13 
where the existing results for the tricritical q — 2 Potts 
model apply 9]. For the 0{n) loop model, the existence 
of tricritical points was further confirmed by transfer- 
matrix analyses for a range of values of n [lol llll |. 
Whereas this work yielded reasonably accurate values for 
some universal parameters, no exact formulas were found 
for these parameters as a function of n. 



Recently Janke and Schakel [l^l reconsidered the nu- 
merical evidence for the tricritical 0(n) model and pos- 
tulated that the conformal classification in terms of the 
Kac formula ^13, 14] of the magnetic exponent, which is 
known to be p = g = m/2 for n = and 1, generalizes to 
other n. 

In order to find the 'missing link', which is the relation 
between n and the conformal anomaly c, a clue is pro- 
vided by the observation that a critical 0{n) model cor- 
responds with a tricritical q — Ti^-state Potts model Q ■ 
It would thus be interesting to bring the tricritical Potts 
model into an even higher critical state. This appears 
to be possible 15, 16] by the simultaneous introduction 
of vacancies and their dual counterparts, four-spin cou- 
plings, into the Potts model. The equivalent loop model 
on the surrounding lattice then appears to have a param- 
eter subspace where the Yang-Baxter equations are solv- 
able ^3 • Out of four branches of solutions parametrized 
by q, one was identified as a branch of tri-tricritical Potts 
transitions [lil fl^ . The conformal anomaly and expo- 
nents of this model are found as a function of q by using 
an alternative representation as a Temperley-Lieb model 

Since the equivalent loop model has loop weight ^/q, it 
seems an appealing possibility that the universal proper- 
ties of the g-state tri-tricritical Potts model match those 
of the tricritical 0(n = y/q) loop model. The confor- 
mal anomaly derived in Ref. 15 is, expressed in n = 
determined by the following equations: 
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yielded scaling dimensions of which 



Xj = (fcj - l)/[2m(m- 
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where we introduced an index j 
is given by Aj — 2 cos[fcj7r/(m -f 
A2 = -1/A, and A3 = -A. 



= 1, 2 or 3, and kj 
1)], with Ai = 1/A, 
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In order to verify the relation with the tri-tricritical 
Potts model, we employ transfer-matrix calculations for 
the loop model on the honeycomb lattice with vacancies. 
These are introduced as face variables located on the el- 
ementary hexagons. They have two possible states: va- 
cant (weight v) or occupied (weight 1 — v). Furthermore 
there is an n-component vector spin Si on each vertex 
i that is surrounded by three occupied hexagons. The 
one-spin weight distribution is isotropic and normalized 
such that J dS — 1 and / dSS ■ S — n. The partition 
function given by jj] thus generalizes to 

Z = ^1,^^(1 -t;)^-^^ / Y[dS,Y[{l + wS,-Sj) (3) 

where the sum is on all configuration variables: £ is a 
subset of the dual lattice and represents the occupied 
faces of the honeycomb lattice. The product over i\£ 
includes all vertices except those of the vacant hexagons, 
iV„ is the number of vacant faces, N is the total number of 
faces, w controls the strength of the spin-spin coupling, 
and (ij) denotes all nearest-neighbor spin pairs. The 
mapping on a loop model proceeds along the same lines 
as m Ref.|3 and leads to the following partition sum 

Z = ^^z;^^(l-^)^-^^w;^™n^' (4) 
c g\c 

where the first sum is on all possible configurations C of 
occupied faces, and the second one over all configurations 
Q of closed loops on the honeycomb lattice that avoid the 
empty faces; is the number of vertices on t/, and Ni 
the number of loops. 

The transfer matrix is constructed for a model 
wrapped on a cylinder, whose axis is parallel to one of 
the lattice edge directions. The unit of length is defined 
as the small diameter of an elementary hexagon. The 
transfer matrix keeps track of the change of the num- 
ber of loops, vacancies, and visited vertices when a new 
layer of sites is added to the cylinder. Its largest eigen- 
value determines finite-size data for the free energy den- 
sity, from which the conformal anomaly can be estimated 
hj. Three more eigenvalues Ai were calculated. These 
determine three correlation lengths and allow finite-size 
estimates Xi(z;, w, L) of the corresponding scaling dimen- 
sions Xi [13. The temperature dimension Xt was esti- 
mated from the second eigenvalue of the transfer matrix, 
and the magnetic dimension Xh from a modified transfer 
matrix with a single loop segment running in the length 
direction of the cylinder. The 'interface' exponent X^ 
was estimated by inserting a column with bond weights 
of the opposite sign. Further details about the transfer- 
matrix technique are given in 

Refs.EHHB 

We parametrize the vicinity of the tricritical fixed point 
by two relevant temperature-like fields ti and ^2, and by 
an irrelevant field u. The associated exponents are yt^ , 
yt2 and i/„ respectively, with yt^ > yt^ . 



TABLE L Tricritical points as determined from the transfer- 
matrix calculations described in the text. The estimated 
numerical uncertainty in the last decimal place is shown in 
parentheses. 



n 


V 


w 


-2.0 


0.3503 (1) 


0.8156 (1) 


-1.75 


0.3649 (1) 


0.8330 (1) 


-1.50 


0.380814 (1) 


0.852082 (1) 


-1.25 


0.3979352 (1) 


0.8726404 (1) 


-1.00 


0.4163568 (3) 


0.8949268 (1) 


-0.75 


0.4360088 (1) 


0.9189617 (2) 


-0.50 


0.4566834 (2) 


0.9446100 (2) 


-0.25 


0.4781475 (2) 


0.9717428 (2) 





1/2 


1 


0.25 


0.521805 (1) 


1.028950 (1) 


0.50 


0.54313 (1) 


1.05812 (1) 


0.75 


0.56361 (2) 


1.08708 (2) 


1.00 


0.5830 (1) 


1.1155 (1) 


1.25 


0.6010 (1) 


1.1429 (1) 


1.50 


0.6175 (1) 


1.1688 (1) 


1.75 


0.6321 (1) 


1.1928 (1) 


2.00 


0.6452 (1) 


1.2145 (1) 



The tricritical point is estimated by simultaneously 
solving the unknowns v and w in the two equations 

X,{v,w,L) ^ X,{v,w,L-l) = X,{v,w,L-2) (5) 

where the functions Xi {i ~ h, t, m) are provided by the 
transfer-matrix algorithm. Expansion of the finite-size- 
scaling function in the vicinity of the tricritical point 
yields that the solution v{L) of Eq. © converges to the 
tricritical value t;(*''') of v as 

i;(L) =„(t") -HaL^"-^'2 +... (6) 

where a is an unknown constant; and w{v, L) similarly 
converges to the tricritical value w^-^"^ . The values Xi{L) 
taken at the solutions of Eq. (|SJ| converge to the tricritical 
scaling dimension Xi as 

X,{L) = Xi + buiy- + ■ ■ ■ (7) 

where b is another unknown constant. We applied this 
procedure both for Xi = Xh and for Xi = Xm to locate 
the tricritical points and to estimate the tricritical expo- 
nents from Eq. |(7|). These calculations were performed 
along the same lines as in Ref. O, but here we use larger 
finite sizes up to L = 14, and moreover we include sev- 
eral values for —2 < n < 0. The results for the tricritical 
points are listed in Table I, together with the estimated 
error margins. The analyses using Xh and Xm generated 
consistent results and allowed us to check the numerical 
uncertainties. 
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TABLE II: Conformal anomaly and tricritical exponents as 
determined from the transfer-matrix calculations described 
in the text. Estimated error margins in the last decimal place 
are given in parentheses. The numerical results are indicated 
by '(num)'. For comparison, we include theoretical values 
obtained from Eqs. Q, and For n < —3/2, the 

temperature exponent Xt becomes complex. 



n 


c (num) 


c (exact) 


Xryi fnum) 


Xm (exact) 


—2.0 


—0 9914 (2) 


—0 9915599 


—0 202 (Vi 


—0 2017990 


-1.75 


-0.9108 (2) 


-0.9109986 


-0.1765 (2) 


-0.1769723 


-1.50 


-0.8196 (2) 


-0.8197365 


-0.15166 (3) 


-0.1516447 


-1.25 


-0.7164 (1) 


-0.7164556 


-0.12596 (3) 


-0.1259301 


-1.00 


-0.6000 (1) 


-6/10 


-0.10001 (2) 


-1/10 


-0.75 


-0.46962 (1) 


-0.4696195 


-0.07410 (1) 


-0.0740955 


-0.50 


-0.32528 (1) 


-0.3252829 


-0.04853 (1) 


-0.0485319 


-0.25 


-0.16799 (1) 


-0.1679953 


-0.023691 (1) 


-0.0236917 

















0.25 


0.17526 (1) 


0.1752630 


0.022111 (1) 


0.0221110 


0.50 


0.35348 (1) 


0.3534792 


0.04224 (1) 


0.0422357 


0.75 


0.52994 (1) 


0.5299489 


0.05999 (1) 


0.0600004 


1.00 


0.70000 (1) 


7/10 


0.0749 (1) 


3/40 


1.25 


0.860 (1) 


0.8589769 


0.0867 (2) 


0.0865052 


1.50 


1.001 (2) 


1 


0.094 (5) 


0.0880192 


1.75 


1.04 (4) 




0.098 (5) 




2.00 


1.05 (2) 




0.10 (1) 




n 


Xt (num) 


Xt (exact) 


Xh (num) 


Xh (exact) 


-2.0 






-0.094 (1) 


-0.0951627 


-1.75 






-0.087 (1) 


-0.0876431 


-1.50 


0.709 (1) 


0.7097847 


-0.0792 (1) 


-0.0790909 


-1.25 


0.4817 (2) 


0.4814739 


-0.0694 (1) 


-0.0693653 


-1.00 


0.4000 (2) 


2/5 


-0.0584 (1) 


-7/120 


-0.75 


0.3445 (2) 


0.3446681 


-0.04593 (3) 


-0.0458895 


-0.50 


0.30390 (2) 


0.3039309 


-0.03199 (1) 


-0.0319828 


-0.25 


0.273220 (1) 


0.2732199 


-0.016645 (1) 


-0.0166435 


0.00 


1/4 


1/4 








0.25 


0.232500 (1) 


0.2324957 


0.017731 (1) 


0.01772952 


0.50 


0.2193 (1) 


0.2192386 


0.03628 (1) 


0.03627658 


0.75 


0.2090 (2) 


0.2088741 


0.05539 (1) 


0.05539746 


1.00 


0.2000 (1) 


1/5 


0.07500 (2) 


3/40 


1.25 


0.193 (1) 


0.1906800 


0.0950 (2) 


0.09549714 


1.50 


0.180 (5) 


0.1684499 


0.12 (1) 


1/8 


1.75 


0.183 (10) 




0.13 (1) 




2.00 


0.184 (10) 




0.15 (2) 





We also obtained finite-size estimates of Xt at the tri- 
critical points thus calculated, and extrapolated these 
data. The results for the conformal anomaly and the 
three exponents are listed in Table II, together with the 
estimated error margins. 

A comparison of the numerical results for the confor- 
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FIG. 1: Conformal anomaly (-|-) and temperature exponent 
(x) of the tricritical 0(n) model vs. n. The data points show 
the numerical data, the curves the theoretical predictions. 
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FIG. 2: Scaling dimensions X^ (x) and Xm (+) of the tricrit- 
ical 0(n) model vs. n. The data points show the numerical 
data, the curves the theoretical predictions. 

mal anomaly with Eq. J^l, as given in Table d strongly 
supports the classification of the tricritical 0(n) model 
as proposed for n < 3/2. Our numerical results for 
Xt agree with X2 in Eq. those for Xm agree with 
Xi. Using the value of the conformal anomaly and m 
as a function n, we confirm that the numerical results 
for the magnetic scaling dimension agree with the entry 
(p = m/2, q = m/2) in the Kac formula 



X. 



p.q 



[p(m + 1) — qmY — 1 
2m{m -f 1) ' 



(8) 



The n > results correspond with branch 1 as defined 
in Ref. [3 and those for n < with branch 2. The 
numerical results and theoretical values of the conformal 
anomaly and the dimensions are shown as function of n 
in Figs, in and 0. 

The expressions for Xi and X2 in Eq. (|2Jl are not re- 
produced by entries in the Kac table, at least not with 
index pairs that are linear in m. This made it difficult to 
conjecture the exact values of X^, and Xt from numerical 
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data alone, even if supported by data for c. This problem 
did not apply to which appears in the Kac table. 

Remarkably, the Potts tri-tricritical branch ends at 
q — 9/4. For g > 9/4 the model is no longer critical and 
the transition probably turns first-order 0, The 
equivalence q — thus yields the result that the tricrit- 
ical 0(n) branch ends at n = 3/2, possibly with a dis- 
continuous transition for n > 3/2. At first sight, the nu- 
merical results for 3/2 < n < 2 may not seem suggestive 
of a discontinuous transition, and allow only a very weak 
discontinuity. But it is clear from Tables I and II that the 
estimated errors tend to increase with n for 77.^3/2, as 
a result of deteriorating finite-size convergence. This is 
consistent with the possibility that an operator becomes 
marginal at n = 3/2, in line with c — 1 (see Table lll|l . 
Therefore, one may expect similar phenomena as for the 
g > 4 Potts model, where the marginal operator leads to 
misleadingly slow finite-size convergence which obscures 
the weak first-order character in a range of q near 4. 

The results presented above apply to the non- 
intersecting loop model. Loop intersections are irrelevant 
in the critical 0(n) model, but they are relevant in the 
low-temperature phase . While the possible relevance 
of such intersections could modify the universal behavior, 
this appears not to be the case for the n = 1 tricritical 
0(n) loop model, since its exponents agree with those of 
the corresponding spin model. 

The scenario sketched above indicates that the criti- 
cal and tricritical 0(n) branches are not connected, and 
does not provide a relation between the tricritical 0(n) 
model and the critical Potts model, such as was recently 
suggested (12i |. 
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